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Based on an explicit verification of the coupling matrix elements be- 
tween the lplh and 2p2h states we propose a new method of selecting the 
most important 2p2h states responsible for fragmentation effects. In this 
way the dimensionality of the problem is reduced, such that the compu- 
tation becomes feasible and the spreading of the strength is realistic, as 
verified by some tests of convergence. Calculations in 208 Pb show that due 
to sizable mixing effects only about 50% of the total isoscalar giant dipole 
resonance (ISGDR) 3fiui strength is located in the energy region between 20 
and 25 MeV. This is the energy region which currently is available in exper- 
iment. Even above 30 MeV we find about 10% of the total strength. This 
indicates that the current experimantal evaluations of the ISGDR centroid 
energy may significantly underestimate its value. 

PACS numbers: 21.60.-n, 24.30. Cz 



In formal terms many interesting nuclear modes of excitation |/), as for 
instance the nuclear giant resonances, are generated by one-body operators 
of the type / = Y, a p fapa\t a P such that 

l/) = />> = £H/|0)l«>, (1) 

n 

The symbol |0) represents the ground state and \n) the spectrum of eigen- 
states in a corresponding subspace. For such excitation modes certain global 
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aspects of the strength function 

S f (E)=Y / Sf(n)5(E-E n ), (2) 

n 

where 

S f (n) = \(n\f\0)\ 2 , (3) 

can thus be described in the subspace of one-particle - one-hole (lplh) 
(|1) = aj,a/j|0}) states generated by the nuclear mean field. In the lplh 
subspace we thus have \n) = J2i c i I n general, however, such states no 
longer remain the exact eigenstates when more complex configurations of 
the npnh-type are taken into account. The lplh components of a much 
larger number of new eigenstates |n) are then spread over many more cor- 
responding new eigenvalues. This is a mechanism of fragmentation. The 
nuclear interaction is predominantly two-body in nature and thus directly 
couples the lplh states to the 2p2h ones only. Therefore, in practical terms 
it is enough [1] to diagonalize the nuclear Hamiltonian 

H = £i4 a i + 4 v ij,kl a i a ] a i a k, ( 4 ) 

i ij,kl 

in the combined space of lplh and 2p2h states. In this equation the first 
term denotes the mean field which in the present work is taken as a local 
Woods-Saxon potential including the Coulomb interaction. The second term 
is the residual interaction with antisymmetrized matrix elements ify- fcj and in 
the following discussion is represented by the density-dependent zero-range 
interaction of Ref . [2] . 

The Hamiltonian matrix then reveals the following structure: 



<l|ff|l'> 


<l|ff|2> 


<1|£|2> 


<2\H\2>> 



Even after such a truncation the dimension of the above matrix is usually 
still much too large to be numerically diagonalized. This in particular holds 
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true for the isoscalar giant dipole resonance (ISGDR) [3]. This 3hu excita- 
tion is typically located in the energy region above 20 MeV and the number 
of relevant 2p2h states is of the order of 10 6 . An interesting conclusion can 
however be drawn from Fig. 1 which shows the distribution of the coupling 
matrix elements (l|ii"|2) between the lplh and 2p2h states in the J w = 1~ 
sector of the 208 Pb nucleus. These states are here generated by the six mean 
field shells (three above and three below the Fermi surface). As it is clearly 
seen from Fig. 1 there is only a very small fraction of the coupling matrix 
elements which significantly differ from zero. This seems to offer an extra 
criterion for selecting the most important 2p2h states. Thus, by setting a 
finite positive threshold value H t h one can select the 2p2h states such that 
|(l|.ff|2)| > H t h for those lplh states which carry the strength. Interest- 
ingly, an explicit numerical verification shows that even a relatively severe 
selection according to this prescription on average preserves a form of the 
distribution of 2p2h states. What in this connection is particularly impor- 
tant is that, as an example in Fig. 2 illustrates, even certain high energy 
2p2h states survive selection and thus a possibility to move the strength to 
higher energies is retained. 

We now evaluate the ISGDR strength distribution using a prescription 
as described above. ISGDR is one of the most interesting nuclear excitation 
modes. This partly originates from the fact that its centroid energy can di- 
rectly be related to the nuclear compression modulus [4]. The corresponding 
one-body isoscalar dipole operator reads: 

/ = r 3 Y! - rrrY u (5) 

where t] = 5(r 2 }/3. The second term in this equation removes the spurious 
center of mass motion component from the operator r 3 Y\ [5]. The resulting 
3Huj strength distribution in 208 Pb on the lplh level is shown in Fig. 3. 
Almost all this strength is located between 20 and 25 MeV. This about 
corresponds to the energy region where the isoscalar dipole strength can be 
identified in the present day experiments on 208 Pb [6]. The picture changes 
however significantly when mixing due to the coupling to 2p2h states is 
allowed. This is illustrated in Fig. 4 which on the three successive panels 
indicates a degree of fragmentation for H t h = 0.4, 0.3 and 0.2 MeV (from 
bottom to top). The number of the corresponding 2p2h states included 
equals 349, 1125 and 4374, respectively. Consistently with our previous in- 
vestigations [7] a specific form of the resulting strength distribution strongly 
depends on many factors and thus also on Hth- However, more global char- 
acteristics, like a percentage of the total strength in certain sufficiently large 
energy windows is much more stable as can be concluded from Table 1 which 
lists such quantities for energies above 25 and 30 MeV, respectively. 

A reasonable convergence of those results, together with a realistic input 
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of the present model, provides quite a convincing indication that one may 
expect about 50% of the total J n = l~ isoscalar 3huj strength in the higher 
energy region, above 25 MeV, i.e., in the region which is dominated by 
many other multipoles and thus this portion of the strength escapes an 
experimental detection. Even above 30 MeV one finds almost 10% of the 
total strength. The present calculations thus suggest that a recent empirical 
estimation [6] of the nuclear incompressibility (Ka = 126±6 MeV) for 208 Pb 
may appear much to low. 

This work was supported in part by Polish KBN Grant No. 2 P03B 140 
10 and by the German-Polish scientific exchange program. 
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FIGURE CAPTIONS 

Fig. 1. Distribution of the coupling matrix elements between the lplh and 
2p2h states for J* = 1" sector in 208 Pb. 

Fig. 2. Energy distribution of all the 2p2h states generated by the six 
major mean field shells (a), (b) corresponds to those 2p2h states (2) which 
fulfil the condition |(l|iT|2}| > 0.3 MeV. The symbol N denotes the number 
of states in the bin of energy equal to 0.2 MeV. 

Fig. 3. Isoscalar 2>%uj dipole strength distribution in 208 Pb calculated in 
the subspace of lplh states. 

Fig. 4. Isoscalar 3Tilo dipole strength distribution in 208 Pb calculated in 
the space of lplh and 2p2h states, for three different values of H t u- 

TABLE CAPTION 

Table. 1. Percentage of the total isoscalar 3huj dipole strength in 208 Pb 
calculated in the space of lplh and lplh states in the energy region above 
25 MeV and above 30 MeV, respectively, for the three different values of 
H t h- The numbers in parenthesis list the corresponding numbers of 2p2h 
states. 



H th =0.2 MeV (4374) 


H th =0.3 MeV (1125) 


H th =0.4 MeV (349) 


S(E > 25 MeV) 


S(E > 30 MeV) 


S(E > 25 MeV) 


S(E > 30 MeV) 


S(E > 25 MeV) 


S(E > 30 MeV) 


46.9 % 


9.2 % 


51.9 % 


8.7 % 


50.3 % 


5.7 % 



